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Abstract. Immortality will remain a fantasy for as long as aging is deter- 
mined by the erroneous biochemical reactions during a particular DNA repli- 
cation. The replication and base excision repair mechanism, associated to 
eukaryotic DNA polymerase-/? enzyme are central to maintaining a healthy 
cell. Here, we give unambiguous theoretical analyses and prove that the ex- 
clusive biochemical reaction involved in a single nucleotide insertion into the 
DNA primer can be efficiently tracked using the renormalized van der Waals 
(vdW) interaction of a stronger type, and the Hcrmansson blue-shifting hy- 
drogen bond effect. We find that there are two biochemical steps involved to 
complete the insertion of a single base (cytosine) into the 3' end of a DNA 
primer. First, the 03' (from a DNA primer) initiates the nucleophilic attack 
on P Q (from an incoming dCTP), in response, 03 a (bonded to P a ) inter- 
acts with H' (bonded to 03'). These interactions are shown to be strongly 
interdependent and require the forming and breaking of P — O and H — O co- 
valent bonds, which in turn imply that we do not need any external energy 
supply. We found that the unique ionization energy values can be used to 
determine why and how a selected set of complex biochemical reactions are 
initiated and executed with negligible errors to maintain healthy eukaryotes. 
Most of the biophysical and biochemical processes, including the biochemical 
reactions addressed here can be readily captured via the generalized vdW- and 
Coulomb-type interactions. 



1. Introduction 

Death comes in many ways and it always does. But on the bright side, it is the 
only biologically viable way to properly put evolution in play. Otherwise, there will 
be no proper speciation due to improper genetic variation, which will defeat the 
purposes of evolution. Rightly so, our attempts here to evaluate the responsible 
biochemical interactions and the reaction pathways leading to the polymerase-/? 
(pol-/3) DNA replication can be related to aging if it is a major contributor. At 
least, our results presented here can provide the essential knowledge required to 
understand pol-/3 and other similar reaction pathways in eukaryotic cells, which 
can be used to "turn-off" genetic related diseases. The above biochemical reaction 
pathways are presently one of the most subtle scientific information needed to 
elucidate what constitute the origin of pol-/3 specificity and fidelity [TJ |2]. We 
chose the eukaryotic enzyme pol-/3 here for two obvious reasons — it is one of the 
simplest and smallest (39 kDa) known DNA polymerase enzyme, and it has been 
well-studied both experimentally and via simulations since the 1970s [3]. 
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Radhakrishnan and Schlick pQ, and Alberts, Wang and Schlick [2] have carried 
out the constrained quantum/molecular mechanical calculations to obtain the rele- 
vant pol-/3 DNA replication pathways with the lowest energies. Here, "constrained" 
means certain variables are computationally made not-to-float, following the crystal 
structure data reported by Batra et al. [4], and from their earlier studies |2]. The 
constrained variables include the position and the numbers of water molecules and 
ions (Mg 2+ , Na + and CI - ), the charge states and the conformational changes of 
the pol-/3 residues (D190 (aspartate 190) and others) and the incoming dCTP (de- 
oxyribocytosine 5'-triphosphate), and so on [2J. The crucial result [TJE] is that the 
initial deprotonation of the terminal DNA primer 5 ~03'H group to a water mole- 
cule occur concurrently with the nucleophilic attack of 5_ 03' on catalyzed by 
Mg 2+ ions, as depicted in Fig. [2j Note that we have used 6+ and 6— instead of 5+ 
and 2— for the phosphorous and oxygen ions, respectively, in order to accommo- 
date the changes of these charges as a result of their changing polarizabilities in the 
presence of changing interaction strengths. In fact, we will implement this notation 
throughout the paper for any element in any molecule without caution. The read- 
ers need not worry about this, because when the time comes, we will identify their 
magnitudes with a set of properly defined inequalities. Any number that carries a 
prime (3' or 5') refers to a specific carbon atom in a given molecule. 

All the computational analyses thus far focused on the "pre- and post-chemistry 
avenues" , which also correspond to the close- and open-conformations of pol-/3 
and DNA complex [H |2j |6] . In these reports [U El El E] , several responsible 
reaction-pathways were proposed entirely based on these avenues, disregarding the 
biochemical reaction route (with the highest specificity and fidelity) in the chemical 
process itself, sandwiched between those two avenues (for more information, please 
refer the Section, Additional notes, given prior to conclusion). Tracking this route 
(within the chemical process) is far more efficient as there can only be a limited 
number of reaction routes to initiate the nucleophilic attack that could insert a 
single base to the growing DNA strand. Therefore, it is technically easier, accurate 
and reliable if we could exploit the information obtained from this chemical process 
to narrow down the processes in the pre- and post-chemistry avenues. Here we make 
use of the ionization energy theory [7] and the energy-level spacing renormalization 
group method [8] to evaluate the van der Waals (vdW) interaction of a stronger 
type, down to the electronic level. We will show that such an evaluation will lead us 
to find the exclusive biochemical reaction (further assisted by the Hermansson blue- 
shifting hydrogen bond [SI [TOl E] effect) with the highest specificity and fidelity, in 
the presence of the correct Watson-Crick hydrogen bonds [T2] . 

2. Theoretical method 

Here, we first introduce the ionization energy theory (IET) and its approxi- 
mation, which are needed to renormalize the vdW interaction, and to determine 
the average atomic energy-level spacing. This is followed by the introduction to 
the renormalized vdW interaction, which will be shown to be responsible for the 
(bio) chemical reactions if this vdW is of a stronger type, much stronger than the 
standard vdW interaction that has a l/R 6 dependence, where R is the separation 
between two nuclei. 

2.1. Ionization energy approximation. Our analyses here use the ionization 
energy approximation, which can be understood in a straightforward manner if we 
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start from the IET-Schrodinger equation [TJ, 



(2.1) ih 



dt 



h 2 

^— V 2 + Viet 
2m 



*(r, t) = H mT *(r, t) = (E ± £)*(r, t), 



where ^(r, t) is the time-dependent many-body wave function, h = h/2n, h denotes 
the Planck constant and m is the mass of electron. The most important parameter 
here is the eigenvalue, Eq ± £, which is the real (true and unique) energy levels for 
a given quantum system (an atom or a molecule or a solid or anything in between). 



Equation (2.1 ) reduces to the standard Schrodinger equation by noting that Eq±£ = 
E, which is exact by definition, and it is not trivial [?J. Here, Eq is the energy levels 
at zero temperature and in the absence of any external disturbances. This means 
that Eq is a constant, while ±£ is known as the ionization energy (or the energy- 
level spacing) where ± refers to electrons and holes, respectively. In this work, — £ 
does not play any role for an obvious reason, but holes (positive charge carriers) 
can play a vital role in biomolecules as proven by Fujitsuka and Majima [13] with 
respect to charge transfer in DNA. There are two formal proofs for the ionization 
energy approximation. The first one is indirect and reads, Hi ET ^(r,t) = (E ± 
C U atter m )*(M) °c (E ± Catom itucnt )*( r ^): available in Ref. [7j, while the second 
formal proof is direct (and is based on logic) and it is associated to the excitation 
probability of electrons and holes within the ionization energy based Fermi-Dirac 
statistics [HI HI)] . 

The ionization energy approximation reads 

icy n\ .-quantum ^constituent 

\ A - A ! Smatter 01 Satom ! 

where Cma'ttci™ is the real energy level spacings of a particular quantum system. In 
the most general sense, £ is defined to be the real energy spacings between energy 
levels such that an electron from one energy level needs to overcome an energy 
cost when that electron tries to occupy another energy level. In molecules, £ refers 
to the energy level spacing between an occupied level (in the highest occupied 
molecular orbital (HOMO)) and an empty level (in the lowest unoccupied molecular 
orbital (LUMO)), which is the molecular energy level spacing. Due to the above 
definition for £, the quantum matter discussed in this work refers to molecules. 
Now, we can readily predict the changes to ^matter™ if we know £atom ltucnt i where 
the values for £atom Ituent can be obtained from the experimental atomic spectra 
available in Refs. [TBI [UJ. For example, £h 2 o oc [2£ h + + £o+]- The ionization 
energy approximation becomes exact if the quantum matter are atoms or ions 
uecaubc ? matter — ? ion . 

In addition, we can also prove why an atomic hydrogen takes up the role of a 
cation in a water molecule, while an oxygen atom ends up as an anion. To prove 
this, we need to verify £ 2 h+ < £o 2 +i which can be confirmed from Refs. [TU H7] 
where £211+ (1312 kJmol -1 ) < £o 2 + (2351 kJmol -1 ). Here, the averaging follows 



Eq. (2.3) (given below). This inequality also proves that the formation of O2H is 
never possible. Therefore, an atomic hydrogen gives up (donates) its electron to an 
oxygen atom, and this process is never the other way round (having oxygen as a 
donor) because there is always a stronger Coulomb attraction between the electron 
from H and the oxygen nucleus. In other words, if £ is larger for atom A compared 
to atom B, then the nucleus of atom A will always attract the valence electrons from 
atom B with a stronger Coulomb force, and not vice versa [7]- This reinforces the 
idea that an atomic hydrogen acts as a cation, whereas, an oxygen atom takes up 
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the role of an anion, in a system consists of H and 0. Consequently, one can use this 
(ionization energy) approximation to evaluate the changing interaction strengths 
between two atoms or ions from different molecules. However, the bonds between 
atoms do not necessarily involve a single valence electron from each atom, and thus, 
we need to average the ionization energies for atoms that donate (or share) more 
than one electron. In this case, the ionization energy averaging follows 



where each subscript j represents one chemical element (Xj) in a particular molecule, 
while i = 1, 2, • • • ,z, where i counts the number of valence electrons coming from 
each chemical element. For example, if X J= i = P such that P forms five bonds with 
its nearest neighbors, then P — > P 5+ , and therefore, 



The first, second, • • • , z th ionization energies (prior to averaging) were obtained 
from Refs. [THIE]- For a water molecule, X J= i = H and X J= 2 = O. It is to be noted 
here that the ionization energy theory was first formally proven to be a renormal- 
ized theory within the energy-level spacing renormalization group method [3 [8] . 
This energy-level spacing renormalization procedure can be related exactly to the 
Shankar renormalization technique within the screened Coulomb potential deriva- 
tion, and his renormalization method is documented here jT5J Q1J1 HD]. Further 
experimental supports for IET and its applications can be found here (2"Tl |2"!Z1 |2"5I 



Moreover, the importance of ionic valence states on the electronic transport 
properties of oxides and magnetic semiconductors was first experimentally proven 
by Mahendiran et al. [3U] , the computational (using the Density Functional Theory) 
proof was given by Mahadevan and Zunger |31) , while the theoretical proofs (using 
IET) have been developed in Ref. [32] . In these reports [SHI [3TJ [33] , the valence 
state of any cation is shown to vary depending on its surrounding. Therefore, 
our averaging procedure explained earlier, properly takes the valence-state effect 
into account. The "surrounding" here refers to the nearest, and the next nearest 
neighbors only, such that these neighbors are directly bonded to this cation. For 
a given biochemical process however, the "surrounding" also needs to consider the 
neighbors that are not directly bonded to a cation or an anion. This additional 
interaction originating from the non-bonding neighbors can be taken into account 
by means of changing interaction strength between two atoms or ions (see below 
for details). 

2.2. Renormalized van der Waals interaction for biochemical reactions. 

Earlier, we have justified why the ionization energy approximation allows us to track 
the interaction strength between atoms from different molecules. Here, we need to 
capture that an increasing attractive strength between two atoms can lead to a 
(bio)chemical reaction. To do that, we need to develop a new vdW formalism that 
allow a stronger Coulomb attraction between two atoms (or ions) for all R > r% + 
where R is the separation between two nuclei, r x and r 2 denote the radii of two 
interacting atoms or ions [33J . This stronger Coulomb attraction has been captured 



(2.3) 




(2.4) 
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using the renormalized vdW formalism, and the relevant formula is [33j 

h 
2 



(2.5) K^siO = V™nb + -hu 







(72 - exp 


[H 



where, Hloq denotes the ground state energy (of a given molecule) in the absence 
of external disturbances, and for T = OK, and V^ Vaals (^) is the stronger version of 
the attractive vdW interaction, after renormalization (33]. Moreover, V^ aals (£) is 
negative guaranteed by the inequality l/y/2— 1 < and the attractive Coulomb 
interaction between an electron and a nucleus is also negative by definition. Here, 
H is the Planck constant divided by 2ir, A = (127reo/e 2 )aB, in which, ob is the Bohr 
radius of atomic hydrogen, e and eo ar e the electron charge and the permittivity of 
free space, respectively. 

The detailed derivation for ^v Vaals (^), its physical insights and why it is suitable 



for both long and short distances are given elsewhere [33]. Equation (2.5) has been 
proven to be stronger than the standard vdW prediction due to the existence of an 
asymmetric polarization between two atoms. Such a polarization exists if the two 
atoms have very different ionization energy values such that ^ tom < £g tom where 
atom A will become a cation, while atom B assumes the role of an anion based 
on IET. Presently [U [6] [3U |35], and usually this asymmetric polarization 
are assumed to be static (by assigning fixed charges to atoms or ions) and there- 
fore, the interaction between any two ions are electrostatic. For example, the red- 
and blue-shifted hydrogen bonds were determined entirely from the step-by-step 
changes to the distance between the supposedly hydrogen-bonded atoms j3H |3S] . 
It is physically alright to assume the charges assigned to the atoms of a particular 
molecule to be static where this assumption can be used to determine the elec- 
trostatic dipole-dipole, dipole-induced dipole and/or the London dispersion forces. 
For example, the above stated interactions can become stronger or weaker entirely 
from the step-by-step changes to the separation between two strongly interacting 
atoms (i?), however, this assumption (by definition) ignores the systematic changes 
to the atomic polarization when R changes. 

The above mentioned atomic polarization and the vdW forces are often treated 
as adjustable parameters in the ab-initio quantum mechanical calculations, and also 
in molecular dynamical simulations |36j . In particular, in any ab-initio quantum 
mechanical calculations, polarization functions are introduced by hand to allow 
some variations to the atomic orbitals in order to evaluate the molecular integrals. 
For example, within the Gaussian basis set formalism, each basis function (indicated 
with a subscript "r" ) represents a normalized linear combination of a few Gaussians, 
namely, the contracted Gaussian [3"T] 

(2.6) \ r ^ djjk\r9ijk 7 

ijk 

where dijk\r is the contraction coefficient, the primitive Gaussian based on the 
Cartesian coordinate (centered on atom b), 

(2.7) g m = Nxly J h z^ exp(-ar^), 
and the Cartesian-Gaussian normalization constant, 



(2.8) N 



2n\ 3/4 r(8a)^'+ fe i!j!fc!" |1/2 
(2i)!(2j)!(2fe)! _ 
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Here, Eqs. (2.6|, (2.7) and (2.8) are designed (not through an ab-initio method) 
to allow polarizations where i, j and k are nonnegative integers, a > is an 
orbital exponent, and Xh, j/t>, £b are Cartesian coordinates such that the nucleus 
b is the origin. When i+j + k = l = 0, one has the s-type Gaussian, while 
the polarizable p- and d-type Gaussians can be obtained for I = 1 and / = 2, 



respectively. Equation (2.6) represents the guessed linear combination of atomic 
orbitals (LCAO), which allow distortions such that an atomic orbital's (AO) shape 
and its center of charge can be shifted upon the formation of a molecule, by requiring 
I > 0. The point is, the above AO with I > are introduced at the outset by 



hand, and these AO-functions (see Eq. (2.6)) are not determined from any ab 



initio quantum mechanical calculation or derivation [371 , and therefore, we can 
now understand why the parameters associated to polarizations were treated as 
variationally adjustable parameters |36j . 

As a consequence, we need to find an alternative way to understand how and 
why the change in the polarization strength can influence the blue- or red shifting 
hydrogen bond (or the generalized vdW interaction), in the presence of another 
ion. These information can be extracted from the ab-initio calculations by im- 



posing some guessed AO-functions (Eq. (2.6)) at the outset, and then treat the 



parameters related to polarizations as variationally adjustable. In our approach, 
the polarizations arc determined directly from IET and its approximation, with- 
out any inbuilt guessed functions or adjustable parameters. On the other hand, if 
one considers an electrostatic interaction, then an induced polarization can change 
the charge magnitude that was originally assigned to an ion, making such an as- 
sumption inadequate to study chemical reactions. Here, we will use the ionization 
energy approximation and the renormalized vdW theory to tackle the polarization 
effect that is responsible to initiate a chemical reaction between two strongly inter- 
acting atoms (or ions) from different molecules. In the following analyses, we will 
show that the changes in the polarization caused by different surrounding molecules 
give rise to the corresponding changes to the partial charges of the reactant atoms. 
These changes in the partial charges assigned to reactant atoms or ions (from differ- 
ent molecules) can lead to a lock-in configuration, which can proceed to a chemical 
reaction, if the configuration has the strongest Coulomb attraction. 

2.3. Benchmark analyses: peptide bond formation (condensation reac- 
tion) . Here, we apply IET and the new vdW formula to evaluate the formation of 
a peptide bond between two a-amino acids, 

(2.9) H 2 NCHRCOOH + H 2 NCHRCOOH ->■ H 2 NCHRCONHCHRCOOH + H 2 0, 

where R denotes an organic substituent, which is irrelevant here. The analyses 
presented here and throughout this paper are rather abstract, but they are suffi- 
cient to achieve our objective to unambiguously determine the correct biochemical 
reaction pathway, without the need to compute any numbers (apart from the ion- 
ization energy averaging). Biochemical processes are indeed complicated, but it 
does not imply that numerical ab-initio quantum mechanical calculations (with 
inbuilt guessed functions and adjustable parameters) are the only way. For exam- 
ple, our approach cannot be used to predict the most probable reaction path and 
its corresponding reaction surface through the calculation of activation energies. 
However, the approach presented here is based on well-defined inequalities derived 
analytically within the ionization energy theory. This method can be readily used 
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to obtain the changes in the partial charges (of the reactant atoms), which will lead 
us to determine the strongest vdW and Coulomb interactions, and subsequently to 
the correct reaction pathway. 

One should note that the issue here is to verify whether the final result, the 
correct chemical reaction, has been correctly predicted or not. However, different 
methods may need different parameters to be calculated and/or predicted along the 
way for consistency checks. These consistency checks are necessary, for example, 
an ab-initio computational method requires the calculation of activation energies 
to find the correct chemical reaction, whereas, IET calls for the evaluation of renor- 
malized polarizations, induced by the renormalized vdW and Coulomb interactions 
to find the correct chemical reaction. This means that, to predict the correct chem- 
ical reaction, IET does not require the calculation of activation energies, whereas 
ab-initio computational methods do not need to determine the strength of interac- 
tions. 

In any case, there is one supreme advantage in using abstract analyses — we can 
extract the precise microscopic mechanisms responsible for a chemical reaction, 
down to the electronic level. For example, knowing the correct chemical reaction 
determined from the calculated activation energies is one thing, but it is entirely 
a different thing to be able to determine what is the microscopic origin or interac- 
tions that initiate a chemical reaction. In other words, why should the interaction 
energy of one particular chemical reaction is the highest (and unique) compared 
to another one? Our approach addresses this question microscopically and ana- 
lytically by means of changing vdW and Coulomb interaction strengths. We now 
proceed to evaluate the chemical reaction responsible for a peptide bond formation 
as benchmark analyses. 

The condensation reaction (given in Eq. (2.9)) is schematically depicted in Fig. [I] 
where the dashed arrows indicate the correct reaction pathways. The alternative 
reaction pathways can involve Ci, On, Hn and Hni (not shown). You are kindly 
advised to refer to Ref. [3T] for a much simpler system, in which the system consists 
of only one water molecule and a surface oxygen atom. Our strategy here is to track 
this condensation re actio n pathway that has the strongest interaction (such that 
I^Waals(OI horn Eq. (2.5) is maximum). This means that, a chemical reaction oc- 
curs continuously until completion or equilibrium (provided the external conditions 
to initiate this reaction are met), guided entirely by the changing instantaneous 
interaction strengths (mostly due to electron-electron and electron-ion types), not 
because the reactants are "consciously aware" that their 

nal destination or products (after the chemical reaction) will have a smaller 
free-energy, or maximum entropy, or energetically more stable. There are five pos- 
sibilities for two a-amino acids to react (ignoring R) , but only one of the reaction 
pathways is responsible for the formation of a peptide bond. We have evaluated 
all the reaction pathways, but we only highlight the correct reaction pathway (in- 
dicated with dashed arrows in Fig. [IJ as a benchmark analysis. Both the correct 
and the incorrect reaction pathways for dipeptide bond formation are similar to 
the reaction pathways discussed for pol-/3 DNA replication (given in the following 
section). 

The polarized charges (<$+ or 8—) for the relevant atoms in an isolated amino acid 
can be determined using Eqs. (2.2) and (2.3), while Table [I] lists all the ionization- 
energy values required for our analyses. From these two inequalities, £n 2 + (2129 
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H 111 



WATER 
DIPEPTIDE I 



Figure 1 . Two amino acids (the carboxylic-group from the amino 
acid 1 and the amine-group from the amino acid 2) react to form a 
dipeptide and a water molecule. We have identified three different 
hydrogen atoms, Hi, Hn and Hni, two different carbon atoms (Ci 
and Cn) and two different oxygen atoms, namely, Oi and On. 
The organic substituent, R is irrelevant in this reaction for obvious 
reasons. Dashed arrows indicate the correct reaction pathways 
(NHi — > COi), while the alternative reaction pathways are between 
these pairs of atoms, OiHni(from amino acid 2) — > OiHni(£rom 
amino acid 1), NHi(from amino acid 2) — > OiHni(from amino acid 
1), OiHni(from amino acid 2) — > COi(from amino acid 1) and 
CHn(from amino acid 2) — > COi(from amino acid 1). The above 
notations imply that the first pair (of atoms) belongs to the acid 
amino 2, whereas the second pair (on the right-hand side of "—>■") is 
from the acid amino 1. These alternatives will not be discussed — 
see text for details. 




kJmol" 1 ) > £ 2 h+(1312 kJmor 1 ) and £ N 4+(4078 kJmor 1 ) > £ C 4+(3571 kJmor 1 ), 
we can conclude that N is negatively charged, N 5_ , Hi is positively charged, Hj + , 
while Ci is slightly negatively charged, Cj - because even though Hn contributes an 
electron to Ci, the negative charge due to this electron (from Hn) is compensated 
by the electron (from Ci) that is polarized toward N to create a covalent bond 
with N. This electron from Ci is actually asymmetrically polarized toward N due 
to £ N 4+(4078 kJmor 1 ) > £ C 4+(3571 kJmor 1 ). Therefore, Hl <5+ > H "£+ and N <5- 
> Cl S— . In addition, a C — C bond gives rise to an isotropic polarization (because 
£c = £c)j an d thus, C does not gain any charge (S+ or 5—) by making bonds 
with another C. However, Cf^~ is positively charged because £c 4 +(3571 kJmol -1 ) 
< £ 4 +(4368 kJmol -1 ), and no charge contribution from Ci to Cn due to isotropic 
polarization. Therefore, we have these inequalities, Cu S+ > Hl <5+ > H "<5+ and 
Oij_ > o u $_^ T ^ c second inequality is due to £ 4H +(1312 kJmol" 1 ) < <^ C 4+(3571 
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kJmol -1 ). Recall here that atoms with a lower averaged ionization energy easily 
gives away its electron, compared to an atom with a larger averaged ionization 
energy. 

Table 1. Averaged atomic ionization energies (£) for neutral and 
charged elements and their respective averaged valence states or- 
dered with increasing atomic number Z . All the experimental ion- 
ization energy values were obtained from Ref. [16]. The notation 
2 nd denotes the 2 nd ionization energy (without averaging), and so 
forth. However, the valence 1+ is also equivalent to the 1 st ion- 
ization energy, while the valence state larger than 1+ has been 



averaged using Eq. (2.3 1 



Elements 


Atomic number 


Valence state / 






Z 


£ number 


(kJmol ) 


H 


1 


1+ 


1312 


-N 


7 


1+ 


1402 


N 


7 


2 nd 


2856 


N 


7 


3 rd 


4578 


N 


7 


4 th 


7475 


N 


7 


2+ 


2129 


N 


7 


3+ 


2946 


N 


7 


4+ 


4078 





8 


1+ 


1314 





8 


2 nd 


3388 





8 


3 rd 


5301 





8 


4 th 


7469 





8 


5 th 


10990 





8 


2+ 


2351 





8 


4+ 


4368 





8 


5+ 


5692 


carbon 


12 


4+ 


3571 


carbon 


12 


5+ 


10423 


P 


15 


1+ 


1012 


P 


15 


2 nd 


1907 


P 


15 


3 rd 


2914 


P 


15 


4 th 


4964 


P 


15 


5 th 


6274 


P 


15 


6 th 


21267 


P 


15 


5+ 


3414 



We can now proceed to evaluate the correct (experimentally proven) reaction 

pathway, NHi(from amino acid 2) interacts with OiHni(from amino acid 1). Hj + 

first interacts with Oj~, which further polarizes the Oj~ electrons (contributed by 

Cjj + and Hjj^) much closer to Oj~. Moreover, the electron contributed by Hj + to 

N 5 ~ is now more polarized closer to N 5 ~ because Hj + now attracts an electron 

from Oj~ originally contributed by Cjj + . These induced polarizations imply that 
Hi-Oi-di£/ + > Oi-Cu£ + and o l -.u 1 -N S t>_ > Hj-n 5 _ ; where Hl • • ■ Oi denotes 
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the vdW attraction within the generalized vdW formalism. Within the hydrogen- 
bond formalism, it is known as the blue-shifting hydrogen bond effect between Hi 
and Oi. Here, the superscript, Cn— Oj represents the covalent bond between Cn 
and Oi. From these inequalities, h i-°i- c " S'+ > °i- c n,5 + an d Oi-Hi-n^//_ > 
Hl ~ N <5— , we can now deduce that ~ interacts more strongly with Cjj + thanks 
to the interaction between Hj + and Oj~. Next we show that the above charge 
inequalities (due to induced polarizations) will lead us to a lock-in configuration, 
which sets the stage for chemical reactions. 

The lock-in configuration is obtained when Hj + interacts with Oj~ and when 
N 5 _ interacts with Cjj + where the strength of these interactions follow Eq. (2.5). 
We can now use the ionization energy approximation to calculate the relevant val- 
ues, namely, £n--c„ cx |£ n + - £c+I = 3f7 kJmor 1 and 6ii-. i k I£h+ _ £o+| = 2 
kJmol -1 , and therefore, the £totai is 317 + 2 = 319 kJmol -1 , which is the highest 
compared to any other alternative combinations of atoms shown in Fig. [I] This 
highest value also gives the highest interaction strength when it is substituted into 
Eq. pJ5| . 

But this is only one-half of the requirement for the condensation reaction to 
proceed. The other half comes from the stronger bonds between Cjj + and Oj~ and 
between N" 5 " - and Hj + . The bond between Cj : + and Oj~ is stronger because the 
electron from Cj : + is polarized closer to Oj~ in the presence of Hj + . Meanwhile, 
the N 5 ~ — Hj + bond has also become stronger because the electron from Hj + 
is now polarized much closer to N 5 ~~ in the presence of Oj _ . In other words, 
these two bonds have become stronger because they are in certain polarized modes, 
which are actually best suited to initiate selective chemical reactions compared to 
unpolarized bonds. For example, it is relatively easy to break a polarized bond 
via an electron and/or a proton transfer. On the contrary, unpolarized bonding 
electrons are more homogeneously distributed in a given molecule, and therefore, 
the charges for atoms within the molecule (5+ and 5—) are either too small or 
negligible to initiate a proper and strong lock-in configuration. In the polarized 
modes however, atoms (in a molecule) have large polarized charges (6+ or 5—), 
which can strongly interact with other atoms or ions from another molecule. 

Note here that thermal energy is not useful to set up a lock-in configuration 
because thermal energy causes non-selective electron polarizations where it will 
excite all the bonding electrons homogeneously in a given molecule. Hence, thermal 
energy can be a useful source of energy required for continuous reactions and/or 
to keep the activation energy lower so that a given reaction (in certain polarized 
modes) can proceed. It is to be noted here that we are not claiming the activation 
energies are (or are not) correlated to bonding energies. We do not require or need 
to use such information. 

We have identified these two bonds, Cjj + — Oj~ and N s ~ — Hj + are stronger 
as a result of selective atomic polarizations. In other words, these bonds are more 
difficult to be broken by applying thermal energy, but they are more easily broken 
if we could lure these charged atoms with opposite charges (with an electron or a 
proton). In particular, the above stronger bonds will lead Cjj + to strongly interact 
with N 5 ~ , while Oj ~~ is strongly interacting with Hj + . All we need now are some 
external energy supply (thermal energy for example) to initiate the simultaneous 
electron and proton transfers, namely an electron-transfer from N s ~ to Cj X + and 
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a proton-transfer, Hj + to Oj ~~ . These transfers are likely to proceed compared to 
other combinations of atoms, which have been predicted from the ionization en- 
ergy approximation presented earlier. This lock-in configuration has the strongest 
attraction, which requires the minimal external energy supply (or the lowest acti- 
vation energy). The simultaneous electron and proton transfers will complete the 
condensation reaction, giving rise to a dipeptide and a water molecule as reaction 
products. 

In summary, the charge (5+ and S— ) redistribution (due to polarization) in 
atoms have given rise to a lock-in configuration, and subsequently have led to 
simultaneous electron and proton transfers in such a way that a new peptide bond 
is formed between N 5 ~ and Cjj + , while the N 5 ~ — Hj + covalent bond is broken. 
On the other hand, a new Oj~ — Hj + bond is formed, while the Cj : + — Oj~ bond is 
broken. We can now move on and use the acquired knowledge to tackle the pol-/3 
DNA replication. 

3. Results and discussion 

Schlick et al. [TJ [5] found the nuclcophilic attack, Al runs simultaneously with 
A2 to B3 proton transfers. In contrast, we show here with unequivocal theoretical 
analyses that a direct deprotonation of H + from 03'H to 03 Q (A3 Q pathway) and 
the Al pathway are activated simultaneously. This concurrent process (with steps 
Al and A3 a ) is chemically driven (initiated by a stronger vdW interaction) with 
the highest specificity and fidelity (see Fig. [2]) . This direct deprotonation is also 
somewhat similar to Abashkin, Erickson and Burt [B], in which, the difference being 
their [6] direct proton transfer is from 03'H to 02 Q , and then to 03 a (after bond- 
rotation) . We use the blue-shifting vdW interaction between electrons from different 
molecules [33;, or from different sections of a bigger molecule, as a determining 
factor for biochemical reactions, which completes one base-pair (CG) to the growing 
DNA strand. The letters, A, G, C and T denote the well-known DNA purine 
(adenine and guanine) and pyrimidine (cytosine and thymine) bases. To avoid any 
confusion, we will always write the name in full to denote the element carbon. 

Of course, the steps, Al and A3„ can only be efficiently initiated in the presence 
of the correct (a) conformation of the incoming dCTP by interacting with the pol-/? 
residues and with the guanine (from the template DNA strand) via the Watson- 
Crick hydrogen bonds, (b) position of Mg 2+ ions and water molecules, and (c) 
position of pol-/3 residues and their interactions with cations, anions and water 
molecules. In fact, these points (a to c) are valid because numerous interactions 
(both quantum and classical types) will come into play in a biologically dynamical 
setting, which have been discussed explicitly by Spyrakis et al. [35] and Jayaram et 
al. [3j5] . Apart from those work, the readers are also referred to an excellent review 
by Joyce [40] for a list of the techniques available to study biochemical reactions 
with respect to DNA polymerases. 

The fundamental information on the generalized vdW interaction leading to the 
insertion of a single base (C) to the growing DNA strand are equally (if not more) 
important (i) to theoretically prove the possibility of replacing phosphorous with 
arsenic in the Californian bacterium, strain GFAJ-1, which was experimentally ob- 
served by Wolfe-Simon and her assemblage of scholars [HI [42] , (ii) to evaluate 
the in-vitro characterization of a DNA nanomachine that maps the changes in the 
spatial and temporal pH, as reported by Krishnan and her colleagues [23] > (iii) to 
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Figure 2. There are two processes involved for a single base (C) 
insertion into the 3' end of the DNA primer, depicted in this 
schematic figure. The first process has been proposed in Refs. QJ[5] 
that consists of five steps, Al to B3. Whereas, the second process 
proven here involves Al and A3 Q occurring simultaneously (follow 
the dashed arrows). The first two major steps (discovered by Rad- 
hakrishnan and Schlick pQ) carry the labels, Al and A2, while the 
other steps (from Bl to B3) are part of the subsequent proton(H + )- 
transfer steps assisted by D190, proposed by Alberts, Wang and 
Schlick 2\ (follow the dotted arrows). The diagonal dashed lines 
on the top right corner denote the possible Watson-Crick hydrogen 
bonds with a guanine (from a template strand) , while the solid di- 
agonal arrow pointing to 03'H in an incoming nucleotide (dCTP) 
will be the new 3' end of the primer (growing) strand, after in- 
sertion. Another diagonal arrow pointing to the left (in the DNA 
primer) extends in the 5' end direction. All the sticks (including 
the isolated sticks inside the ring) and the diagonal solid arrows 
represent the covalent a bonds. The conformations and the posi- 
tions of the DNA primer, D190, dCTP and Mg 2+ are not to scale. 
The number 190 refers to the position of an Aspartate residue at- 
tached to a pol-jS enzyme. 
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study DNA crosslinks [44], and the binding affinity between a ligand and a pro- 
tein |45j . (iv) to understand how a female sex hormone (progesterone) chemically 
activates sperms into "hyperactivity" by increasing the Ca 2+ ions concentration 
(inside sperms) j3H [17] , (v) to expose the selective DNA-dependent protein kinase 
inhibition [48], and the role of elecrosaics in DNA compacion [49], (vi) to examine 
how DNA sequences (■ • • ATCG- ■ • ) and base-pair interactions play the essential role 
in protein recognition [501 1511 152j and in sequencing human genomes |53| , respec- 
tively, (vii) to theoretically evaluate the possibility of exploiting G-quadruplexes 
(a four-stranded guanine structure) as a therapeutic target in oncology to develop 
anticancer drugs [54] . (viii) to study the initiation of human colon tumor that 
grows like a cancer due to pol-/3 DNA replication error [55 , and (ix) to understand 
the mechanism responsible for protein-protein, protein-DNA and other interactions 
between biomolecules [56] , 

3.1. Biochemical reactions in the chemical stage. In this work, we track the 
electron- electron and electron-ion interaction strengths via the changing energy- 
level spacing (presented earlier in the theoretical method section) during a partic- 
ular biochemical reaction. Interestingly, the band structure and energy levels of 
proteins have also been suggested to be useful to understand proteins electronic 
interaction [57]. Here, we will track the above interaction strengths so as to un- 
derstand why and how the nucleophilic attack proceeds with a unique path. In 
fact, there is only one approach that is well-defined and allows one to expose the 
chain of events leading to this exclusive path occurring at the critical point, which 
was first introduced in Ref. [21 to capture the origin of a single water-molecule 
splitting on a MgO surface. The formal proofs that justify this approach are given 
elsewhere [58l 159] . Add to that, these interaction strengths that we will be tracking 
are directly related to the ionization energy theory (IET) and the renormalization 
group method, which are also mathematically well-defined [7J[S]. 

Furthermore, our analyses can be made quantitative by invoking the ionization 
energy approximation, which gives us the license to claim that the molecular energy 
level spacings are directly proportional to the atomic energy level spacings of all 
the atom exist in a given molecule [15, 60 . The atomic ionization energies of all the 
atoms in a molecule is of course different from the molecular ionization energies of 
that particular molecule. They are different because the atomic energy levels will 
have to transform (due to wave function transformation) during a chemical reaction 
to form the molecular energy levels [58j[59]. However, the energy level spacings of a 
molecule containing certain atoms has been proven to be proportional (not equal) 
to the averaged atomic energy level spacings of all the isolated atoms present in 
that particular molecule [7J. 

What we need to do now is to first understand the chain of events giving rise to 
the validity of the steps, Al and A2 within IET. When the negatively charged 5 ~03' 
comes closer to P Q (see Fig. [2|, the s ~03' — H l5+ bond is further reinforced, which 
is the effect we discovered in Ref. [2T] ■ The reason is that the electrons contributed 
by H and carbon to create the 03' — H and 03' — carbon covalent bonds are now 
being strongly bound to 03' in the presence of P„ + , which pulls these electrons 
further away from carbon and H (from now on, we denote this specific hydrogen 
as H'). The s ~03' — 5+ H' covalent bond is necessarily stronger due to induced 
electrostatic interaction between ^"OS' and P Q . The reason is that the bonding 
electron (contributed by H') is now more polarized toward 5_ 03' due to P Q , and 
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this will necessarily result in stronger 5 ~03' — <5+ H' covalent bond. Logical proof: 
the above covalent bond has to get stronger, or weaker, and definitely, its strength 
cannot remain the same in the presence of P^ + . Here, all other conditions from the 
pre-chemistry avenue are relatively weak. They are weak because in the chemical 
stage, these external effects coming from the pre-chemistry avenue do not directly 
involve in the chemical reaction. Of course, the molecules, ions and the residues 
from the pre-chemistry avenue do set the stage for the chemical reaction to proceed, 
but they do not directly involve in the nucleophilic attack as reactants. Now, let 
us assume the said bond ( <5_ 03 / — 5+ H') gets weaker in the presence of P* + , and 
this is only possible if the bonding electron from H' is weakly polarized toward 03'. 
Therefore, one can immediately deduce that this leads to a contradiction because 
the bonding electron contributed by H', has to be more polarized toward 03' due 
to the attractive induced interaction between P* + and s ~03'. 

Consequently, it is energetically inconvenient to break this stronger 03' — H' 
covalent bond by brute force (by thermal energy) , such as the electronic excitation 
of the (03' — H) bonding electron, which is now strongly polarized toward P^ + . 
Therefore, it would be easier by allowing the oxygen (0^^-p) from a nearby water 
molecule to form the required hydrogen bond (between Owat and H') that can 
easily blue shift (via the Hermansson blue-shifting hydrogen bond effect) due to an 
attractive interaction between 03' and P Q . Let us rephrase the last sentence due 
to its extreme importance and also because of its far-reaching effect in biological 
systems. The stronger covalent bond between 03' and H' is due to the attractive 
interaction between 03' and P a , and is actually more easily breakable, not by brute 
force, but when it is subjected to an induced electrostatic interaction provided by 
the 0^ VAT from a water molecule. This means that, the original covalent bond 
between 03' and H' is not easily breakable by a water molecule in the absence of 
the attractive interaction between 03' and P Q . The readers are kindly advised to 
remember this point, or recall the theoretical method section as we will not stress 
it again when used later. 

Now, back to Hermansson, we know that her blue-shifting interaction can initiate 
the Radhakrishnan-Schlick deprotonation (step A2). Unfortunately, due to this 
deprotonation (losing H' to a water molecule), one now requires some external 
energy from an unknown source (this could be some additional heat or interaction) 
to overcome the activation energy and initiate the nucleophilic attack (Al). In other 
words, by having H' bonded to 03', one has the strongest interaction between 
the bonding electron (coming from 03' — H') and the nucleus of P*+ leading to 
smaller activation energy. To see why and how this activation energy is smaller, 
we just need to wait a bit longer because we need to introduce and build the 
arguments properly. Anyway, the other subsequent proton transfers will not be 
captured with IET because these latter paths, Bl to B3 are similar to the hydrogen- 
bonding mechanism stated above, and they are not unique. They are not unique 
because even though these steps, Al, A2 and Bl to B3 require minimal total energy 
compared to other pathways (that involve different residues and so on), one can 
always manipulate any number of pathways to further reduce the minimal total 
energy. For example, we can have a direct deprotonation from the H 2 0- - -H' 
to one of the oxygen in the outgoing pyrophosphate without involving D190, by 
manipulating the "constrained" parameters. Here, the dashed lines, " " means 
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a hydrogen bond, " — " denotes a covalent bond, and "=" is a double covalent bond 
used in text. 

In contrast, we prove here the existence of a unique reaction pathway that con- 
sists of two simultaneous steps, (Al) s ~03' attacks P 5+ while (A3 Q ) 5 ~03 Q coun- 
terattacks + H' (from 03') that has the highest specificity and fidelity. We also 
explain why and how these steps are strategically interdependent that make them 
exclusive. But before we dig deep, let us first properly introduce the ionization 
energy approximation relevant to this work. The origin of these charges, 6— and 
5+ are due to their respective polarization, which can be understood from the av- 
erage atomic ionization energy (£) values |22j . Here, large ionization energy value 
implies large energy needed to polarize or to excite the electron from its initial 
or ground state. For example, one can obtain the average ionization energy val- 
ues from Table[l][Il] and write £ 2H+ (1312 kJmor 1 ) < £ 2+(2351 kJmol" 1 ), and 
therefore, H and O are indeed positively and negatively charged, as they should be. 
This inequality proves that oxygen can never be a cation in water molecules. Using 
the same arguments, we can obtain the inequality, £p < ^o, which is also true for 
all individual (1 st , . . ., 5 th ) ionization energies (see Table u\ to justify the type of 
charges (5+ or 6—) carried by these elements, P* + and O in the 5'-triphosphate 
group. 

Up to this point, it is true that we have no constrained parameters to han- 
dle because we deliberately ignored them, not because ignorance is bliss, but to 
keep all the analyses exact and unambiguous within the chemical process, sand- 
wiched between the pre- and post-chemistry avenues. The complex processes and 
their numerous possibilities involved in these avenues necessitate the introduction 
of constrained parameters. To avoid such parameters, the chemical process is iso- 
lated, which is theoretically valid because the interactions during this closing and 
opening stages, do not directly involve in the biochemical reactions that insert a 
base to the growing DNA strand. Thus, in the following paragraphs, our strat- 
egy is to first implement this isolation correctly, and then discuss the possibility 
of carrying out the pre- and post-chemistry investigations subject to the exclusive 
biochemical reaction pathway. 

Apart from the well-known nuclcophilic attack (Al pathway), there are four 
other possible chemical reactions that correspond to four oxygen atoms bonded to 
P a , namely, P a — 03 a — P, P Q =02 a , P Q — O — carbon, and P a — 0~ (see Fig.[3|V to 
D), in which only one of them should occur simultaneously (with Al). Apparently, 
we need to eliminate any three of them by considering each in turn, and then 
prove A3 a has the highest specificity and fidelity. We have to warn the readers 
that the following analyses are entirely analytic and are inevitably subtle, plagued 
by different magnitudes of polarization for the element, P Q , and the four oxygen 
atoms bonded to P Q in the presence of 03', H', P (valid only for the oxygen in 
Pq, — 03 a — P) and carbon (valid only for the oxygen in P a — O — carbon). 

From the ionization energy approximation (given earlier), we can understand 
the validity of the assigned charges in P 5+ — O s ~ — carbon 5 " 1 ", and consequently the 
steps, Al and A3 a are also correct from the induced Coulomb interaction between 
5+ and S-. Of course, this induced interaction has a deep microscopic origin, which 
will activate the required chemical reactions to extend the DNA strand. But it is 
sufficient to consider only the partial charges now (without going into the details), 
to rule out the three out of the four steps listed above, and to find the crucial 



16 



ANDREW DAS ARULSAMY 



USM carbon 

i_03-'°^° 

° u-<-03a O 



EH 



carbon 
H'<--0 ^02« 

o 



o^\ 



o 



H'<- 
\ 



|-03' 

o 



carbon 

•02a > 

03a O 
\ ^O" 



a 



o 



carbon 
H'<-O s y 

>D"0- 



o 



0' 



o 



FIGURE 3. The above schematic diagrams are sketched to depict 
the nucleophilic-attack region, selected and magnified from Fig. [2] 
for different types of counterattack (counter to the nucleophilic 
attack) between oxygen (bonded to P a ) and H'. The dotted arrows 
indicate the induced electrostatic interactions between 6+ and 8— 
charges for different oxygen atoms bonded to P^ + , for instance, 
(A) s ~03 a (also bonded to another P), (B) 5 ~02 Q , (C) s ~ (also 
bonded to carbon) and (D) 0~. Note here that the nucleophilic 
attack (between s ~03' and P*+) remains the same. The molecular 
conformations and their positions are not to scale. 



oxygen atom (bonded to P Q ) that interact with the strongest induced interaction 
with H'. Once the correct oxygen is identified, we can then go on and exploit the 



existence of the renormalized vdW interaction of a stronger type (Eq. (2.5 1) that is 
responsible for the simultaneous nucleophilic attack Al, and the counterattack A3 a . 
The nucleophilic- and the counter-attack depicted in Fig. [3]A_ to D are due to an 
induced attractive electrostatic interaction between 5_ 03' and P^ + . This induced 
interaction is a prerequisite to form a new covalent bond ( 5 ~03'— P*+), which can 



be obtained from the stronger vdW interaction (see Eq. ( |2.5| )). Using Eq. (2.5), 
we can go on to untangle each chemical step (sandwiched between the pre- and 
post-chemistry avenues) involved during the base, C insertion (from dCTP), and 
find the crucial oxygen (bonded to P Q ) that strongly interact with H' during the 
nucleophilic attack. 

3.2. Al and A3 a steps. The moment s ~03' unleashes its polarizable electron 
on P* + , 5 ~03 a strikes back instantaneously, attacking <5+ H'. As a matter of fact, 
these nucleophilic- and counter-attacks are simultaneous to the extent that one can 
also view it in the reverse attacking order. But this is incorrect. The above original 
order-of-attack (nucleophilic-attack invites the counterattack) can be understood 
by noting Fa S+ > H 5+ and 03 6— > ° 3a 6—, which imply a stronger attractive 
interaction between s ~ 03' and P^ + . The first inequality is obvious with four oxygen 
atoms (forming five covalent bonds) with P Q . Whereas, the latter inequality is due 
to |K(P„ + ) +£(P 5+ )] = 3414 kJmor 1 , 4[£(carbon 4 +) + £(H+)] = 2442 kJmor 1 , 
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and therefore, \ j£(P 5 Q + ) + $(P 5+ )] > 5 [£(carbon 4 +) + £(H+)]. Here, the averaging 
follows Eq. ( |2.3[ ). In other words, the electrons contributed by carbon and H are 
of lower ionization energies, and therefore more strongly polarized toward ''"OS'. 
While the electrons responsible for the partial charge in s ~03 a are coming from 
the less polarizable P (see Fig. |3]A_). Here, P is the least polarizable element (with 
five bonds) compared to carbon (with four bonds) and H (with fixed one bond). 
All bonds are predominantly covalent-type, unless stated otherwise. 

The chemical steps depicted in Fig.[3]A_ are also strategically interdependent, (i) 
The nucleophilic attack will first reinforce the 03' — H' bond, however, this stronger 
bond will be easily breakable only if H' is subjected to s ~03 a . (ii) Conversely, the 
counterattack will only make the P a — 03 a bond stronger, but again, this bond is 
also easily breakable if P* + is allowed to interact with s ~03' . As a consequence, 
we can observe that the two steps, (i) and (ii) nicely correspond to Al and A3q,, 
and they are precisely interdependent. For example, both (i) and (ii) when proceed 
simultaneously, cause these bonds, P a — 03 a and 03' — H' to break, which add C 
to the growing DNA strand due to the formation of P Q — 03' bond. On the other 
hand, a PP^ molecule leaves the scene after the formation and breaking of 03 a — H' 
and P Q — 03 Q bonds, respectively. 

One should not be misled into thinking that the above nucleophilic- and counter- 
attack leading to chemical reactions, are not related to electron transfer after the 
initial electron delocalization. For example, any mechanism involving the breaking 
and forming of bonds will have to first deal with the electron polarization. In 
particular, easily polarizable electron will lead to high probability electron transfers 
or enhanced electron delocalization. Therefore, by knowing the relative magnitudes 
of the electron polarization, one can readily infer the probability of the so-called 
electron-transfer rates and/or the probability of electron delocalization from one 
atom to another. As a consequence, in our analyses, we do not require explicit 
calculations to obtain the electron-transfer rates via the molecular orbital or the 
valence bond theories to determine which bond is likely to be breakable, or which 
one is not. Instead, we can infer the bond forming and/or breaking information 
directly from the atomic polarizability. 

Apart from that, we also stress here that we did not falsely neglect the entropic 
contribution embedded in the Gibbs free energy formalism. On the contrary, we 
do not require to know the Gibbs free energy, enthalpy (reaction energy) nor the 
entropy of reaction in our analyses because as we have stated earlier, we made the 
choice to follow the electrons that strongly interact with other electrons coming 
from any nearest neighbor atom. In other words, we follow the path that likely to 
have the strongest electron-electron and electron-nucleus interactions in accordance 
with Ref. [33| . These strongest interactions will require minimal activation energy 
for a particular reaction to proceed. Again, we are not implying that activation 
energies are (or are not) correlated to bonding energies. 

It is clear that the process consisting of Al and A3 a steps is a type of lock-in 
process due to strong interdependence that can be made exclusive with the highest 
specificity and fidelity, if only we could systematically rule out the other oxygen 
atoms bonded to P Q . This interdependence is actually initiated by a stronger vdW 
attraction between 03' and P Q , and between 03 Q and H' with further assistance 
from the Hermansson blue-shifting hydrogen bond effect. The formation of these 
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bonds, P a — 03' and H' — 03 a are assisted by the Hermansson blue-shifting hy- 
drogen bond effect, at least for the bond, H' — 03 Q . Whereas, the formation of 
P Q — 03' bond is due to the Hermansson-like blue-shifting hydrogen bond effect. 
The readers are referred to Ref. [5T] for some straightforward details on the ap- 
plications of the Hermansson-like bonds. Consequently, these steps do not require 
large external energy supply to break the above-stated bonds. 

The process shown in Fig. [3)3 consists of the Abashkin-Erickson-Burt deproto- 
nation steps [5J. In fact, the 02 a is the only close contender against the A3q, step 
explained above. The other oxygen atoms can be ruled out quite easily. Recall that 
the nucleophilic attack (step (i)) introduced earlier remains the same here, and will 
remain so for all other oxygen atoms (see Fig. [3p and D). The instantaneous coun- 
terattack here is between 02 Q and H', which will strengthen the 02 a —P a double 
bond that should be easily breakable only if is subjected to s ~03' . There is 
no difference in the partial charge between 5 ~02 a and s ~03 a (° 3a 5— = ° 2a 5—) 
because both electrons in each case are contributed by P. However, 02 a =P Q is a 
double bond, hence the single P a — 03 Q bond is in fact twice more easily breakable 
than the stated double bond. Unlike the breaking of P Q — 03 Q bond in the A3 a 
step, the interaction energy from the polarized 03' (during the nucleophilic attack 
given in Fig. [3j3) is distributed equally to both bonds, until one of the double bond 
is broken. In contrast, for the A3 Q step, the same interaction energy is fully avail- 
able to break the P Q — 03 Q bond alone, reinforced by the counterattack between 
03 Q and H'. Thus, the process indicated in Fig. [3^3 is inferior to Al and A3 Q steps 
in terms of the highest specificity and fidelity. 

In Fig. [3p, we require the O bonded to carbon and P Q to counterattack H'. In 
this case, the partial charge in the O atom are contributed by carbon and P, 
whereas two P contribute two of their electrons to s ~03 a , implying °S— > ° 3 "S— 
(recall the averaging presented earlier with carbon contributing a more polarizable 
electron, compared to P). Due to this inequality, the O — carbon bond is the one 
that tend to blue-shift more than the required O — P Q bond. Therefore, the O — P Q 
bond cannot take advantage of the counterattack between s ~0 and |5+ H' when P^ + 
is subjected to s ~03' (because the O — carbon bond will be fortified, instead of the 
O — P Q bond). As such, this third process can also be ruled out as it clearly lacks 
the required specificity and fidelity. 

The logic required to exclude the final oxygen (see Fig. J3p) from the list is 
the easiest of them all. Here, 0~ counterattacks H' without any disturbance to 
the P^ + — 0~ bond. Consequently, we can immediately see that the nucleophilic 
attack cannot proceed until the formation of the P a — 03' bond simply because this 
would require an additional core electron from P a with very large £ (see Table [I]) 
to form the 6 th bond between P^ + and s ~03' . This is definitely unlikely, and so is 
the process suggested in Fig.[3p. 

In summary, it is an unequivocal claim that the steps, Al and A3 Q have the 
highest specificity and fidelity within the biochemical reaction stage sandwiched 
between the pre- and post-chemistry avenues, which has been pointed out much 
earlier. Having said that, the next stage of research can be directed to fine-tune 
the constrained parameters to include the conformational changes to the dCTP in 
such a way that the steps, Al and A3 Q can be allowed to occur. Such steps may 
need some three-dimensional conformational changes to the dCTP, but without 
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disturbing the Watson-Crick hydrogen bonds, where these hydrogen bonds will 
remain intact as one can confirm it from Fig. [2] 

3.3. Additional notes. We have ignored the effects coming from the pre-chemistry 
avenue, for example, due to Mg 2+ ions, pol-/3 residues (D190 and probably others), 
and water molecules. They do contribute directly in the pre- and post-chemistry 
avenues, but they do not act as reactants in the chemical stage. Hence, in the chem- 
ical stage, we can indeed ignore the contribution made by these molecules and ions. 
In fact, we can logically rule out the effects coming from the pre-chemistry avenue — 
any ions or charged atoms (attached to nearby molecules or residues), which are 
not reactants, will affect all the nearby reactant atoms (the correct as well as the 
incorrect ones) equally via the electrostatic interaction. This means that, non- 
reactant molecules do not discriminate (not selective) when they interact via the 
Coulomb interaction (be it repulsive or attractive) with the reactant atoms (directly 
or indirectly) . Interestingly, the effects coming from the pre-chemistry avenue have 
also been ruled out in Ref. [34 using some molecular dynamical simulations. The 
results reported in Refs. [TJ [5] emphasizes the importance of pre-chemistry effects, 
but ignore the contribution of different oxygen atoms bonded to P a in the chemical 
stage because there is no counter-attack in their calculations. This means that, 
their calculations in the chemical stage do not consider the interactions between H' 
and other oxygen atoms bonded to P Q because H' is lost to a water molecule. 

In any case, we have given unambiguous analyses as to why and how the Al and 
A3 Q chemical steps have the highest specificity and fidelity, and are precisely inter- 
dependent (with lock-in steps), which can be made compatible with contributions 
from Mg 2+ ions, pol-/3 residues (D190 and probably others), and water molecules 
in the pre- and post-chemistry avenues. Our approach guarantees that the lock-in 
steps are only possible if H' stays bonded to 03' during the nucleophilic attack. 

As a final remark, this work can also be treated as a proof for the correctness of 
Baskaran hypothesis, "all basic condensed matter phenomena and notions mirror in 
biology" |61j . For example, the same polarization effect discussed herein to explain 
pol-/3 DNA replication, has been applied previously (via the atomic polarizability) 
to prove the existence of carrier- type transitions in oxides |22j . 

4. Conclusions 

In conclusion, we have actually attempted to tell an unambiguous and a self- 
consistent tale of how the electronic structure (with distinct energy-level spacings) 
may influence aging in humans using abstract analyses. These influences are found 
to be indirect where the unique ionization energy values determine why and how 
a selected set of complex biochemical reactions are initiated and executed with 
negligible errors to maintain healthy eukaryotes. Here, we have considered only 
one type of enzyme known as the polymerase-/?, and a specific biochemical pro- 
cess involving the insertion of a single C base (from dCTP) into the growing DNA 
strand. Nonetheless, the above influence (due to ionization energy) is proven to be 
significant and far-reaching, such that the stronger renormalized vdW interaction 
and the Hermansson blue-shifting hydrogen bond effect are applicable to all bio- 
chemical processes. For example, these two interactions are responsible to set up 
the locked-in configuration that has the strongest Coulomb attraction. 

We have proven that there is only one unique biochemical process, consisting 
of two steps with the highest specificity and fidelity. These steps are strongly 
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interdependent to the extent that they can be regarded as lock-in steps, occurring 
simultaneously. In particular, the nucleophilic attack starts when the 03' uses 
its polarizable electron to attract the positively charged P a by means of a stronger 
vdW interaction. This step leads to a stronger 03' — H' bond, which becomes easily 
susceptible in the presence of 03 Q . At the same time, the counterattack between H' 
and 03 Q fortifies the P Q — 03 a bond, which can be easily weakened by subjecting 
P Q to 03', and this is nothing but the nucleophilic attack that is being unleashed as 
we speak. Obviously, these (nucleophilic- and counter-) attacks complement each 
other rather nicely in such a way that they lead to the formation of new 03' — P a 
and H' — 03 a bonds, and the breaking of P Q — 03 Q and 03' — H' bonds. In addition 
to this almost perfect matching, note that the bonds that are being formed are of 
the same types, as the ones that are being broken. This implies a negligible penalty 
in activation energy, or may not even require any external energy supply at all, to 
overcome the activation energy. 
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